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We study the ground-state properties of a dilute gas of strongly-interacting fermions in the 
framework of the coupled-cluster expansion (CCE). We demonstrate that properties such as 
universality, opening of a gap in the excitation spectrum and applicability of s-wave approx- 
imations appear naturally in the CCE approach. In the zero-density limit, we show that the 
ground-state energy density depends on only one parameter which in turn may depend at 
most on the spatial dimensionality of the system. 
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1. Introduction 

One of the major endeavors in modern physics is the quest to isolate the under- 
lying physics in a given system against the background of irrelevant complexity. 
Arguably, this challenge is tied into the ability to investigate systematically the 
predictions of a given theoretical Hamiltonian model. Separating model features 
from artifacts of the theoretical approximations involved, represents the key to the 
quantitative understanding of the physics of strongly- inter acting systems. Unfor- 
tunately, first-principles calculations are presently possible only for small systems. 

To achieve predictive power in many-body theory, it is important to develop a 
theoretical framework that can provide a systematic way of improving over mean- 
field theory and to demonstrate conclusively the numerical convergence of the 
results with the order of the approximation. For strongly-interacting systems of 
particles, traditional perturbation theory fails to converge fast enough for practical 
purposes. Hence, the desired nonperturbative character of the theoretical frame- 
work to be employed. In this context, we submit that the coupled-cluster expansion 
approach to solving the many-body Schrodinger equation represents the best hope 
of extending present state-of-the-art first-principles calculations to the realm of 
systems with large number of particles. 

It is also important to identify problems involving strongly-interacting sys- 
tems of particles that can provide the test bed for many-body theoretical ap- 
proaches and approximations. One such system is the infinite (matter-like) many- 
body system composed of spin-1/2 fermions interacting via a zero-range, infinite 
scattering length contact interaction. The problem of finding the ground-state 
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properties of this system is referred to sometimes as the "George Bertsch prob- 
lem" [l|, [3, 0, 0] and is of particular interest in astrophysics in connection with 
the equation of state for neutron matter. This problem has been revisited re- 
cently with the advent of experimental studies in ultracold fermionic atom gases of 
the crossover from the regime of Bardeen-Schriffer-Cooper (BCS) weakly-bound 
Cooper pairs to the regime of Bose-Einstein condensation (BEC) of diatomic 
molecules 0, 0, 0, S, i, m [HI, E, [3, Q, [1^ • 



The "Bertsch problem" was originally intended as a challenge parameter-free 
model of neutron matter at subnuclear density. In cold atom physics, the regime 



of interest is known as the "unitarity limit" [16|, ll7|, i.e. the limit near Feshbach 
resonances where the s-wave scattering length, oq, of two atoms with different spin 
components is much larger than the inter-particle distance (/cplool >> !)• Here, 
kp denotes the Fermi momentum of the gas, which is conventionally related to the 
total density of particles, pQ, by the noninteracting Fermi gas formula 

where the momentum integral is performed over the interior volume of the Fermi 
sphere, and a denotes the spin component of the fermion, i.e. cr = it ^. 

The "unitarity limit" corresponds to the BCS to BEC crossover in dilute ultracold 
fermionic atom gases and can be reached by modifying either the s-wave scattering 
length or the system density. For an interaction-induced crossover, the unitarity 
limit corresponds to the singularity in the scattering length and the limit is the 
same when approached with positive or negative scattering length. In this limit, 
the correlations are deemed to be significant and the system is expected to exhibit 
universal behavior, independent of the shape of the potential and dependent only 
on the particle density. 

The importance of correlations in the ground state of dilute fermionic matter in 
the unitarity limit is measured by the numerical value of the ratio 

C = e/eo , (2) 

where e and £q denote the ground-state energy densities of the interacting and 
noninteracting systems, respectively. A close upper-bound to the value of ^ was set 
by the quantum Monte Carlo (QMC) study performed by Carlson et al. [18], which 
gave the value ^qmc = 0.44ib0.01. In contrast, the "universal" curve describing the 



BCS to BEC crossover in the standard BCS variational picture [2l|, |22] derived by 
Leggett [13] gives the mean-field numerical value (,mf = 0.59, which suggests that 
beyond-to-leading order effects are responsible for a change of about 25% when 
compared to the mean-field solution. Recent theoretical and experimental values 
for ^ are summarized in Ref. [19] . 
In this paper, we will argue that by applying the coupled-cluster expansion 



(CCE) [2d, |2J] formalism to the Bertsch problem, we can demonstrate conclu- 
sively that the resulting CCE equations are consistent with expectations based on 
purely heuristic arguments, such as universality and the possibility of a gap in the 
single-particle (s.p.) excitation spectrum. The latter is important for the ability 
to capture the superfluid properties of the system, if present. We will show that 
in the combined unitarity and zero-density limit, the ratio ^ depends on only one 
parameter which in turn may depend only on the spatial dimensionality of the 
system. 
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2. Coupled-cluster formalism 

We will review now the basics of the CCE approach: Consider the many-body 
Schrodinger equation, H = E \^), with a two-body Hamiltonian density oper- 
ator, H = T V. 

In second quantization, operators are expressed in terms of the creation and 
annihilation operators in momentum representation, c^(q) and c(q), subject to 
the canonical relations: 

{ct(qi),c(q2)} =<5(qi-q2), {ct(qi), ct(q2)} = {c(qi), c(q2)} = , (3) 

where spin (and isospin) degrees of freedom are implied. The coordinate and mo- 
mentum space representations of the creation operators are related via the sym- 
metric Fourier transform, 

V;(x) = / c(q) e''^ '' . (4) 



(2^)3/2 

Next, we introduced the particle and hole operators 

a(q)=c(q)0(g-A:F), b(q) = ct(q) ^(feF - g) , (5) 
such that the physical vacuum obeys the relations 

a(q) m=0, b(q) 1$) = . (6) 
With these definitions, we have 

ct(q) = at(q) + b(q) , c(q) = a(q) + bt(q) . (7) 

In order to derive the CCE equations, we begin with the following ansatz for the 
many-body wave function, 1^*) = |^>), where l^") is the physical vacuum, and S 
is the many-body cluster correlation operator defined as 

S = Si + S2 + S3 + • • • , (8) 
where S^ gives rise to the npn/i-configuration contributions in |^), i.e. 

Sn =A / d^pi . . . [ d^p^ [ d^K ■ ■ ■ f d^h 



>A;f J>kF J <kw J <k 



X 5'„(pi,- ••ki)aj,^---at^bj^^---bj^^ . (9) 

By construction, the many-body wave function, l^'), obeys the normalization con- 
dition, (^1$) = 1. 

In the case of an infinite system, the physical vacuum, |$), is represented by the 
noninteracting Fermi gas system. Because of translational invariance arguments, 
the s.p. representation of the vacuum is introduced in terms of plane-wave wave 
functions and the amplitudes Sn{pi, • • • p^; k„, • • • ki) in S^, see Eq. ([9]), satisfy the 
property 

n n 

Sn{Pl, • • • Pn; kn, • • • ki) OC 5^ Pi - ^i) ■ (10) 

i=l i=l 
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It follows immediately that the Ipl/i-correlations term in S vanishes, i.e. Si = 0. 

Using the CCE ansatz, the Schrodinger equation can be written in normal or- 
dered form, as 

l[e-W]J$)=e|$), (11) 

where the subscript c indicates the creation part of a normal-ordered operator. For 
c = 0, Eq. (jlip gives the expression for the ground-state energy density, which, 
in the case of a two-body Hamiltonian, i.e. if V = V^^^, includes at most con- 
tributions due to 2p2/i-correlations. For c 7^ Eq. pT]) gives rise to a system 
of nonlinear equations that must be solved self-consistently for the amplitudes 
Snip 1) ■ ■ ■ Pni ■ ■ ■ ki) ill S^. In the case of an infinite system, the equation corre- 
sponding to c = 1 is identically zero because Si = 0. Furthermore, approximations 
to the CCE system of nonlinear equations are based on the idea that the Sn am- 
plitudes are small as long as the relative distance between particles is much larger 
than their average distance. Hence, approximations such as S^ = for n > N, 
are equivalent to truncating an expansion in powers of density [25]. Therefore, in 
the dilute limit, we can disregard contributions due to S^ with n > 3, and Eq. ([8]) 
gives simply S = S2. We conclude that, for an infinite matter system of fermions in 
the dilute limit, we only need to solve the equations for c = (the energy-density 
equation) and c = 2 (the equation for the amplitudes in S2). 



3. Interaction model 

Without loss of generality, we consider here the same two-body potential used by 
Carlson et al. [18|] in their QMC study of the unitarity limit, i.e. 

1 

V{r) = -- . (12) 

m cosh (ar) 

By construction, the reduced wave function at zero energy in this potential is 

ito(^) = ripoir) = tanh(ar) . (13) 

It follows that the potential (|12p corresponds to an infinite s-wave scattering length, 
oo = — 00. We find the range of the potential (fT2]) is 



ro= / drMr) - Mr)] = -, (14) 
Jo " 

where we have introduced the notation, (pQ = I — UQ^r), to denote the difference 
between mo('') and its asymptotic behavior. Finally, we note that the dilute limit 
is achieved by taking the limit ar^ — > 00. Here, denotes the unit radius defined 
as 

d-TT 

Poyrf = l. (15) 



We obtain = l/i-ykp), with 7 = [4/(97r)] 
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Figure 1. Diagrammatic representation of the energy of an infinite two-body Hamiltonian system. Here, 
vertical single and double lines depict hole and particle states, respectively, i.e. states with momenta 
k < k-p and p > A:fi respectively; filled circles indicate the kinetic energy operators, whereas horizontal 
wiggles and thick lines indicate the two-body potential operator and the 2p2h-correlation amplitudes, 
S2 (pi , P2 ; k2 , ki ) , respectively. 



We will write the energy density of fermionic matter in terms of the ratio 

V{y/k^) _ 107^ {avs? ^^g^ 
Eo 3 cosh^ [y 7(ars)] ' 

where we have introduced the notation y = kpr. With our choice of potential, see 
Eq. (112p . the ratio ^ depends on the product arg, where a is a measure of the range 
of the interaction and is related to the Fermi gas density, as described above. 
We note that Carlson et al. have carried out QMC calculations for ar^ values of 
12 and 24, and have noted that the changes in ^ between these two values of ar^ 
were negligible within statistical errors. 



4. Energy density in the Bertsch problem 

We recall that, in the case of an infinite matter system, we have Si = 0. Then, the 
energy density reads 

£ = — (^.le^^He^l^) = — (^'||t + V+ rV,S2l||$) . (17) 
Po Po ^ ^ 

This equation is exact and is illustrated diagrammatically in Fig. [TJ 



4.1. Vacuum expectation of the kinetic energy 

The first term in Eq. (I17p represents the expectation value of the kinetic-energy 
operator in the vacuum state. As stated above, in the case of infinite matter, the 
vacuum is the noninter acting Fermi gas. Hence, we have 




which gives the familiar result, Eq = (3/5)eF; with ep = kp/{2m). This allows one 
to write the exact expression for the ratio ^, as 

e = l + ^(^l|v+[V,S2]||$>. (19) 

Therefore, the quantity (^-1) measures the departure from the noninteracting Fermi 
gas result. This is given by the sum of two contributions: i) the expectation value 
of the two-body interaction operator, V, in the vacuum state, and ii) the expec- 
tation value of the commutator [V, S2] in the vacuum state. (We note that, by 
construction, we have (<I>| [V, S2] = (<I>| V S2|<I').) 
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a) 

Figure 2. (Color online) The integration domain of the center-of-mass coordinate, K, depends on the value 
of the magnitude of the relative coordinate, k, as follows: a) < fcp — fe, b) fcp — fc < < + k'^ , 

c) Jkl + A;2 <\K. The shaded areas indicate the Fermi spheres (A; < fcp). 



4.2. Vacuum expectation of the potential 

We discuss now the contribution of the two-body potential, V, alone. We have 

[2S + 1) ■ 

5=0,1 



(20) 



where S denotes the total spin of a pair of fermions, k and K represent the relative 
and center-of-mass coordinates of the pair, respectively, and the anti- symmetric 
matrix element of the interaction is defined as 

(q'lV^^-'^^lq) = (q'|V(^)|q) + (-l)^(-q'| V^^^ |q) . (21) 

For a central potential, V = V{r), the matrix element of the interaction reads 

(q'|V(^)|q) = J dV V{r) e'('i-'i'>'" . (22) 

In Eq. (120p . Dk indicates the integration domain for the variable K, which is 
illustrated in Fig. [2l Using 



IV. (2vr)3 

where we have introduced the notation, 



(23) 



(24) 



we obtain 

($|V|$) = 4 / dk k^B{k) / drr^ [l^(r-)/eo] [2- jo(2A;r)] . (25) 

Po^o Jo Jo 

Here j(,{z) denotes the regular spherical Bessel function of rank ^. 

We note that, using the power expansion of jo{z), we can show that the radial 
integral in Eq. (j25p can be calculated only in terms of the even momenta of the 
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potential, (r^"F(r)), i.e. 

foo 2^.2 

/ drr^ V{r) [2 - jo{2kr)\ = {V{r)) + — (rV(r)) , (26) 

Jo 3 

with 

{r^^Vir))= [drr^]r^''V{r). (27) 
Jo 

We evaluate the ratio V{r)/eo using Eq. (jl6p . Then, the density dependence of 
the expectation value of the interaction in the vacuum, see Eq. (j25p . is illustrated 
in Fig. [3l We notice that in the zero-density limit, this contribution to the energy 
density vanishes. Using the explicit form of the interaction, Eq. (jl6p . we can show 
that the momenta of the potential, (r^"y(r)), decrease in size with increasing n, at 
fixed density. So, according to Eq. ([26]). in the dilute limit the contribution due to 
($1 V|$) depends only on {V{r)), which is independent of the shape of the potential. 
Therefore, independent of the explicit choice of the inter-particle interaction, we 
obtain that the zero-density limit of the ratio ^ in the CCE approach is obtained 
as 

C^eo = l + ^(^l[V,S2]|^). (28) 
Poeo 

Hence, the universal character of this result in the unitarity limit. Moreover, the 
quantity (Co-1) represents a direct measure of the importance of correlations in the 
ground state. 

4.3. Correlations effects 

We turn now to the calculation of the correlations effects. We begin by recalling 
that according to Eq. (jlOp we can write 

52(pi, P2; k2, ki) = Z(p, k; K) S\P - K) . (29) 

Here, k is confined to the volume of the Fermi sphere, i.e. k < kp, p spans the entire 
space and the center-of-mass momentum K is confined to the domain illustrated 

^<(D|(V/£o)|4>> 




Figure 3. Vacuum expectation value of the two-body potential, with the inter-particle density given by 
the infinite s-wave scattering length potential introduced by Carlson et al. [T8| . Here, the dilute limit is 
recovered for ars S> 1- 
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Figure 4. (Color online) Diagrammatic representation of the S2 equation in the case of a dilute infinite 
matter system of spin-1/2 fermions. The four terms boxed on the first line read as ^2 (pi , P2; k2, ki) (ep^ + 
ep2 ~ ^fei ~ ^fe2)' where '^p(k) a-re the renormalized s.p. energies given in Fig. |5] 



in Fig. [2j Then, we obtain 

($|[V,S2]|«') 



(30) 



2 ^ 

5=0,1 



/■d3pZ(^'-)(p,k;K) (k|V(^)|p). 



Next, we perform the partial wave decomposition of the amphtude Z('^)(p, k; K), 
i.e. 



Z(5)(p,k;K) = J2 zfjAP^k-K) { \v^'^\p)®Y^'^\k)\ QY^^\k)\ , (31) 

where we denote by Y^^\q) the spherical harmonic of rank £ that depends on 
the angular coordinates of the vector q. To calculate the expectation value in 
Eq. (j30p . we need to solve the S2 equation illustrated diagrammatically in Fig. SI 
This equation is nonlinear, as expected because of the nonperturbative character 
of the CCE approach. For illustrative purposes, we note the terms boxed in on the 
first two lines in Fig. H] lead to the familiar approximation 



'S'2'^'''''(Pi,P2;k2,ki; 



(piP2|V(^'")|k2,ki) 

£pi + £p2 ^fcl €fc2 



+ 



(32) 



where ep{fc) are the renormalized s.p. energies illustrated in Fig. [5j If instead one 
uses the noninteracting s.p. energies, then the S'2(pi, P2; k2, ki) amplitudes are 
singular when all momenta are located on the Fermi sphere. Fortunately, as seen 
from Fig. [5l the renormalized s.p. energies include renormalizations for interaction 
(terms 2 and 3) and correlations (term 3) effects. We note that terms 1 and 2 have 
the same signs for both particle and hole states, whereas terms 3 have opposite 
signs and give rise to a gap in the s.p. energy spectrum. 

In order to reveal the true nature of the correlations effects in (j28p , it is necessary 
to introduce the 2-hole function, defined as \2t 



(33) 



xS(r) = / d'p e'P" ^(^)(p,k;K) 



This is equivalent to Fourier transforming the p degrees of freedom of the ampli- 
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tudes Z(p,k;K). The associated partial-wave-expansion components of the 2-hole 
function obey the relation 

poo 

xfX^ir, k-K)= dpp'hM) zf^L^V. k; K) , (34) 
Jo 

and have the anti-symmetric correspondent 

x5:l(r, k;K) = {l + l [(-)^+^^ + i-f^'^] } x[%,ir, k; K) . (35) 
We obtain 

^ ($1 [V, S2] 1^) = 4^ E (-1)' ^2lTT (36) 



POEO TT 



x/ dkk''B{k) I drr'[V{r)/eo]h{kr)xf,i{kr;0), 
Jo Jo 

where the spin component 5 obeys the selection rule £ + S = even. From Eq. 
Xe^e^^ii^i k', K) obeys an equation of the form 



m 



= -(-l)^V4vr(2^i + 1) [1 + (-1)^+'^] V{r) j,^{kr) 5,,,, 5lo , 

which is similar to the radial part of a nonlinear scattering problem. 

In particular, we are interested in the 2-hole function partial-wave component, 
Xuo^^ ikr;0), that enters Eq. ([36j) . In order to write the general form of xTm\kf) 
it is useful to limit first the discussion to the case of the s-wave approximation 
of i.e. £ = = 0, and consider the power expansion of Xoo'o'*(^^' 0)' 

xS^(^;0) = E^"^"' (38) 

where z = kr. In Fig. [6] we depict the density dependence of the contributions 
to the vacuum expectation value, ($| [V, S2] We notice that powers with 
n < — 2 lead to divergent contributions, so only n > — 1 are allowed in Eq. (|38ll . 
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Figure 5. (Color online) Diagrammatic representation of the renormalizcd s.p. spectrum in the case of a 
dilute infinite matter system of spin-1/2 fermions. The first line gives the pariicie-state s.p. energies, Cp, 
whereas the second line gives the hole-state s.p. energies, efe. In either case, the renormalized s.p. energies 
include the noninteracting s.p. energies, ep = p'^/{2m) and ej. = k^/2m., corrected for interaction (terms 2 
and 3) and correlations (term 3) effects. We note that terms 1 and 2 have the same signs for both particle 
and hole states, whereas terms 3 have opposite signs and give rise to a gap in the s.p. energy spectrum. 
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Figure 6. (Color online) Contributions to the vacuum expectation value, (•I'l [V, S2] I"!"), corresponding to 



powers z 



-2, ■ ■ ■ 2, in the power expansion of X^oqq\^)- Here, we have 



: kr. 



Unlike the solution of a typical scattering problem, the solution Xoo'o^(^^' 
allowed an irregular part at the origin for fixed k. We also find that the powers 
corresponding to n > vanish in the zero-density limit. Therefore the value of 
(^0 — 1) from Eq. (f28]l is proportional only to the coefficient c_i. 

In general, the 2-hole function partial- wave component, X'^iQ^kr;^)), that enters 
Eq. (I36p has the form 



XuQ\kr; 0) = Ai,ji,{kr) + B^neikr) Sio + Xm"'{kr; 0) . 



(5,a), 



(39) 



Solving the S2 equation is equivalent to finding the expansion coefficients Ai and 
Bi and the corrections Xuo^\t^''')- Similar to our discussion surrounding Fig. [H we 
can show by direct computation that the only allowed irregular contribution to 
($1 [V, S2] 1$) corresponds to ^ = 0. Hence, the irregular spherical component in 
Eq. (IMl). 

In Fig. [7] we illustrate the density dependence of the contributions to the vac- 
uum expectation value, ($| [V, S2] |$), corresponding to the no(fcr) and the ji{kr) 
components, with ^=0,1,2. We note that all je{kr) components vanish in the zero- 
density limit, whereas the only nonzero contribution corresponds to the irregular 
s-wave spherical component, nQ{kr). At low but finite density, e.g. avg > 5, the 
contributions due the je{kr) components decrease with increasing values of i. For 
arg > 10, the dominant contributions correspond to the s-wave approximation, 
i.e. £ = 0. Therefore, the value of (^0 — 1) from Eq. ([28l) is proportional to the 
coefficient Bq. In this limit, the only scale present in the problem is the number of 
spatial of dimensions. Thus, the coefficient Bq may depend at most on the spatial 
dimensionality of the system. 



5. Conclusions 

To summarize, in this paper we report results of a formal study of the ground-state 
properties of dilute fermionic matter in the unitarity limit, carried out using the 
CCE framework. In this approach, we are able to demonstrate properties of the 
strongly-interacting fermionic matter such as universality, the presence of a gap in 
the excitation spectrum and applicability of s-wave approximations in the dilute 
limit. We note that these results were obtained in an ab initio fashion and did 
not invoke explicit gap parameters to be optimized variationally such as it is done 
in the BCS mean-field picture. This is particularly important, because it assures 
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- 0.20 L 

Figure 7. (Color online) Contributions to the vacuum expectation value, (<1?| [V, S2] I*!"), corresponding to 
the regular and irregular spherical Bessel components for (£=0,1,2). Here, we have z = kr. 

that a CCE-based description of the density-induced BCS to BEC crossover in 
the unitarity hmit is indeed possible. In the zero-density hmit, the ground-state 
CCE equations show that the ground-state energy density depends on only one 
parameter, which in turn may depend only on the spatial dimensionality of the 
system. In three spatial dimensions, the departure from the Fermi-gas energy- 
density result, measured by the numerical value of (^^0-1)5 is proportional to the 
expansion coefficient of no(A:r), the only allowed irregular piece of the hole function, 

xff^^\kr). The /10/e-function formalism introduced in Ref. ^] is the key ingredient 
that allows one to obtain this result. 
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